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This paper presents a method for estimating time delay margin for model-reference adaptive control of 
systems with almost linear structured uncertainty. The bounded linear stability analysis method seeks to 
represent the conventional model-reference adaptive law by a locally bounded linear approximation within a 
small time window using the comparison lemma. The locally bounded linear approximation of the combined 
adaptive system is cast in a form of an input-time-delay differential equation over a small time window. The 
time delay margin of this system represents a local stability measure and is computed analytically by a matrix 
measure method, which provides a simple analytical technique for estimating an upper bound of time delay 
margin. Based on simulation results for a scalar model-reference adaptive control system, both the bounded 
linear stability method and the matrix measure method are seen to provide a reasonably accurate and yet not 
too conservative time delay margin estimation. 


I. Introduction 

Adaptive control is a potentially promising technology that can improve performance and stability of a conven- 
tional fixed-gain controller. The ability to accommodate system uncertainties and to improve fault tolerance of a 
control system is a major selling point of adaptive control since traditional gain-scheduling or fixed-gain control meth- 
ods are viewed as being less capable of handling off-nominal operating conditions. Nonetheless, these traditional 
control methods tend to be robust to disturbances and unmodeled dynamics when operated as intended. 

In spite of the advances made in the field of adaptive control, there are several challenges related to the imple- 
mentation of adaptive control technology in safety-critical systems. The absence of the verification and validation 
methods of adaptive control systems remain a major hurdle to the implementation of adaptive control in safety-critical 
systems . 1 - 2 This hurdle can be traced to the lack of performance and stability metrics for adaptive control which poses 
a major challenge that prevents adaptive control from being implemented in safety critical systems. The development 
of verifiable metrics for adaptive control will be important in order to mature adaptive control technology for use 
in operational safety-critical systems. Of these, stability metrics of adaptive control are an important consideration 
for assessing system robustness to unmodeled dynamics and exogenous disturbances. In one aspect of verification 
and validation, a control system is usually certified by demonstrating that it meets an acceptable set of requirements or 
specifications for stability margins, among other things. Herein lies a major challenge for verification and validation as 
there is no existing standard tool for stability margin analysis of nonlinear adaptive control. The lack of stability met- 
rics for adaptive control is viewed as a technology barrier to developing certifiable adaptive control for safety-critical 
systems . 1 - 2 
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Classical LTI control systems are certified by demonstrating that they meet specifications for stability margins 
among other things. Typically, certification requirements for flight control systems, such as MIL-F-9490D, are often 
addressed in terms of phase and gain margins. These margins are used for LTI control laws to provide robustness 
or safety margins in a control system design as a safeguard against unmodeded effects and unstructured uncertainty. 
While the gain margin concept has been extended to adaptive control, 2 3 the phase margin concept does not easily lend 
itself to adaptive systems due to the inherent nonlinearity in adaptive control. Strictly speaking, phase margin for 
adaptive control in a global context as in the LTI framework is not possible. However, in a local context, it may be 
possible to consider phase margin in an approximate local sense, keeping in mind that even a standard gain-scheduling 
control system design may also contain nonlinear effects due to gain scheduling. 

Phase and gain margins can be used as stability metrics in adaptive systems under some circumstances. One 
possible use would be when an adaptive control process is terminated by switching off the adaptation, essentially 
freezing the adaptive parameters, or when the adaptive signal converges to a steady state value. Some methods of 
approximate phase and gain margin analysis for adaptive control have been proposed that could be used without 
turning off the adaptation. One method for analyzing stability margins is based on a LMI approach by transforming 
the nonlinear adaptive control into a linear parameter varying form. 4 In another approach, it is proposed to define a 
LTI system that bounds the closed-loop adaptive system and then evaluate the phase and gain margins for the bounded 
LTI system in a local time window. 5-6 A potential benefit of this approach is that the adaptation can be “driven” or 
adjusted on-line to meet an approximate phase margin specification to improve the time delay margin of the closed- 
loop system. 7 Both approaches in 4 and 5-6 use similar system error dynamics. 

Time delay margin has been viewed as a more readily accepted metric for relative stability of nonlinear control. 
While time delay margin is a suitable stability metric for adaptive control, a current challenge is that there is no well- 
established analytical tool for computing the time delay margin. Other methods for estimating the time delay margin 
have been proposed. One such method applies a Pade approximation to approximate a time-delay system. 8 The Pade 
approximation transforms the original time-delay system into a higher order system without the time delay that can be 
analyzed by the Lyapunov method to estimate the time delay margin. However, the Lyapunov method with the Pade 
approximation yields highly conservative estimates of time delay margin even for a simple scalar adaptive control 
system. The discrepancy between the time delay margin estimated by this method and the numerical evidence from 
simulations is at least three orders of magnitude. 8 

Despite the fact that new theoretical methods are being developed for computing time delay margin, they are still 
not ready to be used in a unified framework like the classical phase and gain margins. The most direct way to compute 
time delay margin is by simulations. The time delay margin is estimated by introducing a time delay at the input of an 
adaptive control system and then adjusting it until the closed-loop system is on the verge of instability. However, for 
adaptive control to be accepted in the future, simulation-based time delay margin computation is not considered to be 
sufficient as long as there is a lack of analytical tools for the same. 

This paper presents a new method for estimating time delay margin for model-reference adaptive control of systems 
with almost linear structured uncertainty. Bounded linear stability analysis method has recently been introduced to 
represent the conventional model-reference adaptive law by a locally bounded linear approximation within a small 
time window using the comparison lemma. The locally bounded linear approximation of the combined adaptive 
system is cast in a form of an input-time-delay system over a small time window. The time delay margin of this 
system represents a local stability measure and is computed analytically by a matrix measure method, which provides 
a simple analytical technique for estimating an upper bound of time delay margin. Based on simulation results for 
a scalar model-reference adaptive control system, both the bounded linear stability method and the matrix measure 
method are seen to provide a reasonably accurate and yet not too conservative time delay margin estimation. 

II. Introduction 

III. Time Delay Margin for Linear Time Invariant Control 

A. Time Delay Margin for a Simple Scalar System 

Consider a scalar time-delay system 

x (t) = ax (t) + bu (t) (1) 

where x(t) : [0,°°) ->R, u(t ) : [0,<=°) -> R, and b > 0. 

The system has a feedback control 

u(t) = —kx(t) (2) 
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with a closed-loop pole s = a — bk < 0. 

The closed-loop system is designed to be robust to a time delay at the input such that 

x (t) = ax (t) + bu (t — t d ) (3) 

where t d is defined as a time delay margin for the original system in Eq. (1). 

The Laplace transform of the closed-loop time-delay system is expressed as 

(s — a + bke~' dS ) x (s) = x (0) (4) 

To calculate the time delay margin t d , a number of approximate methods can be used. Consider the following 

1. Taylor’s Series Approximation: 

The term e~‘ dS can be expanded using the Taylor’s series as 

e -‘ dS = 1 _ t d s + ] -tjs 2 - ... ( 5 ) 


Expressing in time domain, one gets 

u (t t d ) — u (t) - t d u (t) + ^ t d ii (0 - • • • (6) 

The first-order Taylor’s series approximation in effect is a finite-difference approximation of a time derivative 

cinr*p 

^ ( 7 , 

td 


Then, the time-delay system can be approximated as 

bt 2 

x (t) = ax ( t ) + bu ( t ) — bt d u (t) + (t ) — • ■ ■ 


( 8 ) 


So, the effect of time delay shows up as time derivatives of the controller. This linear equivalent system is 
now conditionally stable. To see this, consider only the first-order approximation of the closed-loop time-delay 
system as 


x (t) = {a — bk)x{t) + bkt d x (t) 


(9) 


which can also be written as 


x{t) 


a — bk 
1 — bkt d 


■x(t) 


( 10 ) 


One can see that the system can be not guaranteed to be stable even if a — bk < 0 since there is an additional 
requirement 

1 - bkt d > 0 (11) 

that must be fulfilled. 

Thus the time delay margin estimate of the system is given by 


1 

bk 


( 12 ) 


This indicates that the time delay must be kept low if the feedback gain is large for the closed-loop system to be 
stable. The result is independent of the parameter a. The Taylor’s series approximation of the time delay term 
does not yield a proper transfer function since the number of zeros is greater than the number of poles. 


2. Pade Approximation: 


Pade approximation is frequently used to approximate time delay 
the form 


P(s) 


ew 

R(s) 


effects by a rational polynomial function of 

( 13 ) 
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where the degree of the polynomial Q (s) is less than or equal to that of R (5) for a proper transfer function 
representation of the time delay. 

Consider the following first-order Pade approximation 


e -t<is = 1 
2 + 


(14) 


The system characteristic equation is 

tjs 2 + (2 — atd~ bktd)s + 2bk — 2a = 0 (15) 

which results in a time delay margin estimate of 


The result now is dependent on all system parameters, 
the order of the Pade approximation. 

3. Lyapunov-Krasovskii Method: 

Stability of time-delay systems can be analyzed using the Lyapunov-Krasovskii method. 910 Consider the fol- 
lowing Lyapunov-Krasovskii functional 

1 

V(x(t)) = x 2 (?)-| — / a' 2 (t)c/t>0 (17) 

td Jt-tj 

The time derivative of V (x(t)) along the solution trajectory is evaluated as 

V (x(t)) = 2 ax 2 (?) — 2 bkx(t)x(t —td) + —x 2 ( ? ) — —x 2 ( t — td ) (18) 

td td 

By completing the squares, one obtains 

V (x(t)) — ^2 a + bk + x 2 (t) + ^ bk — — ^ x 2 (t — td) — bk [x(f) +x(t — td)] 2 (19) 


a + bk 

The accuracy of the estimation increases with increasing 


Since bk > 0, the time-delay system is uniformly stable if the following inequalities are satisfied 


2,& H - bk H — 0 

td 

bk— — <0 
td 

The solution of the inequalities is feasible if a < 0 and bk < —a. This yields 


1 

2 a + bk 


<t d < 


bk 


( 20 ) 

( 21 ) 


( 22 ) 


The result of the time delay margin estimate based on the Lyapunov-Krasovskii method is generally non-unique 
and is dependent upon the selection of the Lyapunov-Krasovskii functional. For example, suppose the following 
Lyapunov-Krasovskii functional is selected 


V(x(t)) 


1 r 

px 2 (f) H — / px 2 (t) dx > 0 

td Jt-tA 


(23) 


where p > 0, then 

V (x(t)) = 2 pax 2 (?) — 2 pbkx 2 (?) +2 pbkx 2 (?) — 2 pbkx (t)x(t — td) + —x 2 (?) — — x 2 (t — td) (24) 

td td 
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which becomes 


V (x(t)) = ( 2 pa + p 2 b 2 k 2 + — ) x 2 (?) + ( 1 — — ) x l it — td) — \pbkx(t) +x(t — tj)} 


P \ .. 2 , 


The solution of the inequalities 


td J \ td 


2pa + p 2 b 2 k 2 + - < 0 
td 

i-^ <o 

td 


yields 


td <P < 


b 2 k 2 


1 + A/l- 


b 2 k 2 


for a < 0 and bk < —a. 


(25) 

(26) 
(27) 


(28) 


4. Lvapunov-Razumikhin Method: 

The Lyapunov-Razumikhin method can be considered as a subset of the more general Lyapunov-Krasovskii 
functional approach. 10 However, a nice aspect of the approach is that it utilizes functions as opposed to func- 
tionals as the main ingredient. The Lyapunov-Razumikhin theorem states that the system is asymptotically 
stable if there exists 77 > 1 and P = P T > 0 such that 

V(x(t))<-e\\x(t)\\ 2 (29) 


where e > 0, whenever 


V (x(t + 0)) < T]V (x(t)) 


for all 0 £ [—tj, 0). 

Consider the following Lyapunov candidate function 

V(x{t))=x 2 (t) 

Differentiating V (x (t)) along the solution trajectory of x(t) yields 

V (x(t)) = 2 ax 2 (t) — 2 bkx(t)x(t — td) 


(30) 


(31) 

(32) 


Recall from fundamental theorem of calculus that 


x(t-t d ) =x(t) 


x{x)di 


n-td 


Then 


(33) 


V (t) = 2cuc (t) — 2bkx 2 (t) — 2bkx (t) / [ax(x) — bkx(x — t d )]dT 

Jt-t d 


<2 (a — bk)x 2 ( t ) +2 bk |x(t)| 


/ [ax(x) -bkx(x — tdj] dx 

Jr-t d 


< 2(a — bk)x 2 (t) + 2\a\bk\x(t)\ f \x(x)\dx + 2b 2 k 2 \x(t)\ f \x{x-t d )\dx (34) 

Jt-td Jt-td 


Since t — t d < X <t, then the Lyapunov-Razumikhin theorem gives 

\x{x-t d ) | < |x(t)| < |x(f)| 


(35) 


Thus 


V(x(t))<2(a — bk)x 1 (t)+2\a\bk\x(t)\ f \x(t)\dx + 2b 2 k 2 \x(t)\ f \x(t)\dx 

Jt-td Jt-td 

= 2 (a — bk)x 2 (t) + 2t d \a\ bkx 2 (?) + 2 tdb 2 k 2 x 2 (?) 

= 2 [a — bk + td(\a\bk + b 2 k 2 )]x 2 (t) (36) 
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So the time-delay system is asymptotically stable if 

a — bk + td[\a\bk + b 2 k 2 ) <0 
Hence, the time delay margin can be found as 


bk — a 
bk(bk + \a\) 


(37) 

(38) 


It is noted that the time delay margin of the time-delay system of Eq. (3) can actually be found exactly by 
computing the system pole s = a + jco as follows: 

<J + jco — a + bke^ td<y e~-’ m ‘ l = 0 (39) 

The system is neutrally stable for (7 = 0 so that the following equations result 

— a + bkcoscotj = 0 (40) 

co — bk sin cotj = 0 (41) 

The solutions of these equations yield the /'co-axis cross-over frequency and the time delay margin as 



The solution actually tends to the Taylor’s series approximation for bk a. Also, there exists a relationship 
between a and bk such that the system is stable, independent of time delay. This occurs when a < 0 and bk < — a . 

Example: Given a = 1 and bk = 2, the Taylor’s series approximation yields t*, = 0.5 sec, whereas the Pade approx- 
imation yields t*, = 0.667 sec. The exact value is t*. = 0.604 sec. Thus, the Pade approximation gives a better estimate 
than the Taylor’s series approximation, but also over-estimates the time delay margin. The Lyapunov-Krasovskii 
method has no solution since a > 0. The Lyapunov-Razumikhin method gives t*, = 0.167 sec. 

Given a = — 1 and bk = 2, the exact value is now t*j = 1.209 sec. The Pade approximation over-estimates the 
time delay margin with ft = 2 sec. The Taylor’s series approximation yields the same estimate of ft = 0.5 sec which 
is independent of a. The Lyapunov-Krasovskii method also provides no solution since bk > —a. The Lyapunov- 
Razumikhin method gives f*j = 0.5 sec, which is the same as the Taylor’s series approximation. In fact, for a < 0, both 
the Taylor’s series approximation and the Lyapunov-Razumikhin method produce the same result. 

Given a = — 1 and bk= \ , the system is stable, independent of time delay. The Taylor’s series approximation yields 
tj = 2 sec and the Pade approximation yields t*, = —4 sec. The time delay margin corresponding to the Lyapunov- 
Krasovskii functional (17) is G = 2 sec and that for the Lyapunov-Krasovskii functional (23) is t d = 7.464 sec. The 
Lyapunov-Razumikhin method gives t^ = 2 sec, which again is the same as the Taylor’s series approximation. 

In these examples, both the Lyapunov-Krasovskii and Lyapunov-Razumikhin methods are quite conservative even 
for a simple linear time invariant scalar system. Relaxation of the conservatism in the Lyapunov-Karsovskii and 
Lyapunov-Razumikhin is possible and usually requires parameter and functional optimization. 

B. Time-Delay Margin by Matrix Measure Method 

Lor a vector time-delay system 

x (t ) = Ax (t) — BKx (t — td) (44) 

where x(t ) : [0,°°) — > R" and A (A — BK) £ C~, i.e., A — BK is Hurwitz, the time delay margin can be found from the 
following characteristic equation 

det (j col - A + BKe~ imd ) = 0 (45) 

The bounds on ft) and t,j can be estimated by a matrix measure method. 9 Defining id as an eigenvalue of a symmetric 
part of a complex matrix such that 

Bi{C) = A, ^ (46) 
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where C € C is a complex matrix and C* is its complex conjugate transpose, then li has the following properties 


fJ-i (C) L M 


( 47 ) 


m = 

max A ; 

\<i<n ' 

fc+c*\ 

OJ 

II 

V 2 ) 

±L(C) = 

min A i 

\<i<n ' 

(c+c*\ 

| = lim 
1 e^O 

K 2 ) 


||/ + gC7||-l 
£ 

l-||/-eC7|| 

£ 


P(C) = -m(-C) 
JI(jC) = -ji_(jC) 
l±(C) <ReA,-(C) <ju(C) 
ImA (C) < ju (— jC ) 


M(C)<||C|| 


jU (C + £)) < ju (C) + /i (£)) 
ju (C + D) > /i (C) + /i (£)) 
JI(C-D)>JI(C)-JI(D) 


H(C-D)<n(C)-n(D) 


(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 


The matrix measure jl affords a simple way to estimate the bounds on the system poles for a MIMO system. 
Lemma 1: The time-delay system (44) is asymptotically stable if the following inequalities hold 


t d < — cos 
CO 


fi(A)+n(jBK) 

\\BK\\ 


co<n(-jA) + \\BK\ 


where ||.|| = ||.|| 9 is the J^-norm. 

Proof: The real parts of the system poles are bounded from above by 
a = ReA, (A- BKe~ ja>td ) < JI (A) +Ji(—BKe~ jmd ) <JI(A)+JI(-BK) |cosffl? rf +~p(jBK) sin cot d \ 
Let 0 < cot d < j, then the time-delay system is stable if cr < 0 which implies 


(59) 

(60) 


(61) 


/i (A) < — /x (— BK ) cos cot d — jl ( jBK ) sin cot d = jl (BK ) cos cot d — /x (jBK) sin cot d 

< JI (BK) cos cotd — /I (jB K) sin cot d (62) 


Upon some algebra, this can be expressed as 

[jl 2 (BK) + Jr (jBK)] cos 2 cot d - 2/T (A ) ju (BK) cos cot d +Jr (A) -Jr (jBK) > 0 
The solution yields a bound on time delay margin t d as 


1 _i 

td < — cos 
CO 


J[(A)Ji (BK) + JI (jBK) \JjI 2 (BK)+JI 2 (jBK) — JI 2 (A) 


JI 2 (BK)+JI z (jBK) 




But 


So 


JI 2 (BK) < JI 2 (BK) +JI 2 (jBK) < HBA'II 


td < — cos 
CO 


1 JT (A) IIBJ5TII +JI(jBK) y/ \\BK\\ 2 -JI 2 (A) i JI(A) + JI(jBK) 


\\BK\\ 


< — cos 
CO 


\\BK\\ 


(63) 

(64) 

(65) 

( 66 ) 
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The imaginary parts of the system poles are bounded from above by 


co = IrnAj (-jA + jBKe~ i<otd ) < p (- jA ) + p (jBKe~ jmd ) 


< B ( ~jA ) + jl (jBK) |cos (Ot d \ + jl ( BK ) sinffl? d (67) 


which can be expressed as 


» <n~JA) + Jp 2 (BK)+J[ 2 (jBK) < Jl (—jA) + ||**|| 


where for added conservatism the less than or equal sign is replaced by the less than sign. 

Corollary: The time-delay system (44) is asymptotically stable independent of time delay if 

jr(A)<||fltf||<-/x(A) (69) 

Proof: The time-delay system is stable, independent of time delay, if 

jl(A)+ jl (—BK) cos cot d + jl (jBK) sin cot d < jl (A)— Jl (BK) cos cot d + jl (jBK) sin cot d 

< M (A) + yV 2 (BK) +jl_ 2 (jBK) = ft (A) + ^p 2 (BK) +JL 2 (jBK) <jl(A) + ||2MT|| < 0 (70) 

This implies 

At(A)<-M (71) 

Note that this condition is in addition to the requirement that A — BK is Hurwitz, which can easily be shown by a 
similar argument that 

P(A) < ||Bjq (72) 


Therefore 


jl (A) < < - jl (A) 


Example: Given 


The bounds on co and t d are computed as follows: 


0 

1 

, BK = 

0 

0 

-1 

1 


0 

2 


jl(A) = l 

P (jBK) = 0 
\\BK\\=2 

co <p (—jA) + llBA'II = 3 rad/sec 

1 _,p(A) K n „„„ 

t d < — cos 1 = — = 0.349 sec 

co BK 9 


The exact results can be determined from 

jco -1 

1 jco -1+2 (cos cot d - j sin cot d ) 


= — co 2 +2 co sin cot d + 1 — jco (l — 2cos cot d ) = 0 


V3 + Vl 

co = = 2.189 rad/sec 

2 


* 2 n 

t d = — —= = 0.478 sec 

3 (y/3 + y/l) 

Thus, the time delay margin estimated by the matrix measure method is reasonably conservative but not overly 
conservative that renders it impractical for design and analysis purposes. 
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Example: Given 


A = 


0 

_ 3 
2 


1 

_ 3 
2 


BK = 


0 0 
1 1 


This time-delay system is stable, independent of time delay, since 


p (A) =0.041 

p (A) = -1.541 
|[BK|| = 1.414 
p(A) < ||5K|| < -p {A) 


C. Time-Delay Margin by Lyapunov-Karsovskii Method 


Stability of time-delay differential equations based on the Lyapunov-Karsovskii method have been studied exhaus- 
tively by many authors. 9, 10 As shown above, different Lyapunov-Karsovskii functionals lead to different results. 
Invariably, the negative-definiteness of the time derivative of a Lyapunov-Karsovskii functional results in a linear ma- 
trix inequality that can be solved for a time delay margin. While this study does not focus on the Lyapunov-Karsovskii 
method, it is instructive to illustrate this technique for estimating time delay margin. 

Lor the same time-delay system (44) with an assumption A (A) £ C , then consider the following Lyapunov- 
Karsovskii functional 

1 

V ( t ) = v T (t)Px (t ) -I — / x T (t) P. v (t) d T > 0 (74) 

td Jt-tj 


where P = P T > 0. 

Evaluating V (t) yields 


V (t) = i ,T (. t)Px(t ) +x T ( t)Px(t ) H — x T (t)Px(t) — — x T ( t — t d )Px(t — t d ) 

td t d 

= x T ( t)A T Px(t ) — x T (t — td)K T B r Px(t) +a ,t (t)PAx(t) — x T (t)PBKx(t — t d ) 

+ — x T (t) Px (: t ) — — x T (t — td) Px (i t — td ) (75) 
td td 


Lor stability, V (t) < 0 is satisfied by the following LMI 


A J P + PA + ±P -PBK 

td 

-K J B J P -J-P 

td 


<0 


(76) 


IV. Time Delay Margin for Adaptive Control 

A. Model Reference Adaptive Control 

Given a nonlinear plant 

x(t) =Ax(t) +B[u ( t ) +f(x(t))+A(x(t),u(t) ,z(t) ,f)] (77) 

z{t)=g(z{t),x(t),u(t),t) (78) 

where x(t) : [0,°°) — > M" is a measurable or observable state vector, z(t) : [0, °=>) — > M. q is an unobservable state vector, 
u(t) : [0,<=°) — > W is a control vector, A G R" x " and B £ W xp are known such that the pair (A,B) is controllable, 

f(x(t)) : R" — > R p is a matched structured uncertainty, A (x(t),u(t) ,z(t) ,t) i s a matched unstructured uncertainty or 

unmodeled dynamics, and g (z(t) ,x(t) ,u (t) ,t) : [0,°°) x R" x R p x R 9 — » R' y represents the unmodeled dynamics. 

The unstructured uncertainty A (x(t) ,u (t) ,z(t) ,t) can be due to numerous unmodeled effects and are generally 
cannot be captured in the control model due to modeling difficulty. In flight vehicles, these effects can represent various 
complex modes of interactions including aeroservoelastic modes, pilot interactions, nonlinear unsteady aerodynamics 
near stall and post-stall, atmospheric disturbances such as sharp-edged wind gusts and wake vortices; just to name 
a few. The tendency for the unstructured uncertainty A (x(t) ,u(t) ,z(t) ,t) is to destabilize a control system since it 
cannot be accounted for in a control design. In a typical control design framework, stability margins are built into a 
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control system to accommodate for the unstructured uncertainty A (x(t) ,u (?) ,z(t) ,t) while the structured uncertainty 
f(x(t)) is kept as small as possible by increasing the accuracy of the control model. The certification for stability of a 
control system is based on meeting well-accepted specifications for stability margins, such as MIL-F-9040D standards 
used for certifying flight control systems. 

The matched structured uncertainty / (x) has a form of 

/(x)=©* T <I>(x) (79) 


where 0* £ W xp is an unknown constant weight matrix that represents a parametric uncertainty, and <t> (x) : R" — > R m 
is a vector of known functions 

The objective is to design a controller that enables the plant to follow a reference model 


Xm (t) — A m x m (f ) T B m r ( t ) 

1 is Hurwitz and known, B,„ £ R nxp is also known, and r (f) : [0,°°) 


where A m £ R" x " i 
with r £ Jzfoo. 

Defining the tracking error as e ( t ) = x m (t ) — x (t ), then the controller u ( t ) is specified by 

u(t) = K x x {t) + K, r (f ) — u a d (x (f )) 


(80) 

£ .Sfoo is a command vector 


(81) 


where K x £ R pxn and K r £ R pxp are known nominal gain matrices, and u a d (x(f)) » R^ is a direct adaptive signal. 
Then, the tracking error equation becomes 


e (f) = x m (t) -x(t)= A m e (t) + ( A m - A - BK x )x(t) + (B m - BK r ) r(t) + B I u ad (x (t )) - 0* T <f) (x (t )) 


(82) 


We choose the gain matrices K x and K, to satisfy the model matching conditions A +BK X = A,„ and BK r = B m . 
The adaptive signal u ad is an estimator of the parametric uncertainty in the plant such that 

u ad {x) = ® T <S>{x) (83) 

where © £ R" txp is an estimate of the parametric uncertainty 0*. 

Let © = © — ©* be an estimation error of the parametric uncertainty. Then the tracking error equation can be 
expressed as 

e(t)=A m e(t) + B@ J (t)<t>{x(t)) (84) 

The system can be designed to follow the reference model with a direct model reference adaptive control update 

law 

©(f) = ~T<t>(x(t))e J (t)PB (85) 

where P = P T > 0 solves the Lyapunov equation 

PA, n +AlP=~Q (86) 


where Q = Q r > 0. 

If A (x (t ) , u {t ) ,z(t) ,t) = 0, then the adaptive law (85) can be shown to be stable and results in e (t ) — > 0 as t — > °°. 
As the adaptive gain T increases, the tracking error further decreases. The upper limit of the adaptive gain is set by 
the sampling frequency of the discrete-time implementation of the adaptive law. Concomitant with the increase in the 
adaptive law is an increase in the high frequency content in the adaptive signal. The estimate © (f ) is essentially a 
nonlinear integral gain since © (t ) can be expressed as 

©(/) = -/’ T^{x{t))e J (t)PBdx (87) 

Jo 

For a LTI system, as an integral gain increases, the closed-loop poles move away from the real axis along the 
imaginary axis. Hence, the frequency of the closed-loop system increases. The phase margin for an LTI system 
generally decreases with increasing the integral gain. Even though &(t) is a nonlinear integral gain, the behavior is 
similar to that of an LTI system. 
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In real systems, because the unstructured uncertainty A (x,u,z,t) is not always zero, the adaptive law (85) is not 
robust and cannot guarantee boundedness of tracking error. Many robust modification schemes can be incorporated 
into the adaptive law (85) to improve its robustness to unmodeled dynamics such as <7 -modification 1 1 

©(f) = -r r<l>(x(f))e T (f)Pfi + c7©(f) (88) 

e-modification 12 

0 (t) = — r <P{x(t))e T {t)PB + n e T (t)PB ©(f)! (89) 

or optimal control modification which has recently been introduced 13 

©(f) = -r <l>(x(f))e T {t)PB — v < t > (x(f)) < J >T (x(f)) © (t)B T PA~ l B (90) 

where c 7,/i, v >0 are tuning parameters in the modification schemes. 

All these modifications invariably introduce new parameters to adjust the adaptive law, for example, cr is such a 
parameter in the above adaptive law. The modification parameters effectively add a damping mechanism to the adaptive 
law to ensure that the adaptive signal is bounded. However, in general there exists a trade-off between performance 
and robustness. In the adaptive law above, increasing the adaptive gain T results in a better tracking performance 
but poorer robustness, while increasing the parameter cr improves robustness but results in a poorer tracking ability. 
Thus, a current challenge in adaptive control design is to be able to select appropriately the tuning parameters that 
can achieve stability and performance specifications. Currently, there is no well-accepted stability and performance 
metrics for adaptive control design and analysis. 

B. Bounded Linear Stability Method 

Global stability is the ultimate requirement for any control system including linear time invariant and nonlinear adap- 
tive control. Global stability is difficult to prove even by the Lyapunov method since it requires a complete detail 
information of a system. Because of uncertainty, it is more tractable to design a control system to satisfy specified sta- 
bility margins rather than attempting to prove that the system is globally stable. Linear time-invariant control systems 
have been designed with the classical phase and gain margins which are well accepted and understood. 

In dealing with nonlinear adaptive control, time delay margin has been proposed by numerous authors as a substi- 
tution for the classical phase margin to indicate relative stability of a nonlinear adaptive control system. Unfortunately, 
there is no existing nonlinear theory for global stability of time-delay nonlinear adaptive control. Attempts to over- 
come this theoretical gap including methods that approximate a time-delay system with a Pade approximation and 
then invokes the Lyapunov method to obtain an estimate of the time delay margin. The global stability requirement 
based on the Lyapunov method inherently results in highly conservative estimates of a time delay margin. To relax the 
global stability, linearization of the nonlinear adaptive control system can provide some local stability estimates using 
the classical linear stability margins. To obtain an equivalent LTI system, the adaptive law can be linearized at a certain 
point in time when the weights have reached their steady state values, usually long after initial transients have settled 
down. However, transient responses during adaptation can be important and can compromise system robustness. 

As an alternative to linearization, the bounded linear stability analysis method has recently been introduced to 
approximate a nonlinear adaptive law with its bounded linear version over a short analysis time window using a 
comparison lemma. The bounded linear approximation of the adaptive law then allows the LTI stability concept to 
be analyzed for the adaptive law locally within the analysis window. The resulting time delay margin estimated by 
the bounded linear stability method using the phase margin and the gain cross-over frequency has been shown to 
have a reasonable agreement with simulation results. In this study, the bounded linear stability analysis provides an 
approximate locally bounded linear version of a nonlinear adaptive law. Using this locally bounded linear system, a 
time delay margin can be estimated using the matrix measure method as derived previously. 

The bounded linear stability analysis method is based on a version of the comparison lemma 14 stated as follows: 

Lemma 2: The equilibrium state y(t) = 0 of the differential equation 

y(0 = -4> T (0rT>(f)y(0 (91) 

where y(t) : [0, °°) — > R, <f> (t ) £ \ : [0, °°) — > R" is a piecewise continuous and bounded function, and F > 0 £ R" x ", 
is uniformly asymptotically stable, if there exists a constant y > 0 such that 

4> T (T)r4>(T)t/T> y (92) 
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( 93 ) 


which implies that y (t) is locally bounded by the solution of a linear differential equation 


i{t) = -yz{t) 


fort G [tj,t, + Tq), where f,- = t,-i + Tq and i = 1,2, . . . ,n — > 

Proof: Choose a Lyapunov candidate function and evaluate its time derivative 

V{t) = \y 2 (t) (94) 

V (t) = -4> t (t) T4> (t )y 2 (t) = — 2<t> T (t) TO (t) V (t) (95) 

Then, there exists y > 0 for which V is uniformly asymptotically stable since 


V(t + T 0 ) 


V (f)exp 




4> : (T)r«5(T)£?T 


<V{t)e~ 2YT ° 


(96) 


This implies that 


Thus, the equilibrium y ( t ) 


exp ^—2^ " <t> T {T)r<t>(x)dT^J <e~ 2Yr ° 

0 is uniformly asymptotically stable if 

1 ft+To 

jrJ <P T {T)T<P{T)dT> y 


provided <J>(f) G ££2 is bounded. 

Then y(t) G ££2 H ££00 since 


v(t- 


K”) — 




Jo 


(97) 


(98) 


(99) 


It follows that 

V(t) < -2yV(t) =>y(t)?(t) < - yy 2 (t ) (100) 

which implies that the solution of Eq. (91) is bounded from above if y (t) > 0 and from below if y ( t ) < 0 by the local 
solution of 

z(t) = —yz{t) (101) 

for? G [ti,ti + To), where to = 0, /,■ = 1 + To, and i = 1,2, .. . ,n — > «>. 

Equation (101) also applies for <I> = <1> (3; (t)) since the condition ( t> (y (t)) G £££ is identically satisfied given that 
y ( t ) G ££2 n££L. This is shown by evaluating V ( t ) as 


dv 

v(t) = y(t)- r = 

dy 

-$ T (y (f )) ro (y (t))y (t ) ^ = -24> T (y (t )) DJ. (y (f )) V (0 

ay 

(102) 

Thus 

dV 0 dy 

V(t) y(t) 

(103) 

Suppose there exists y such that 

XI a. 

IA 

(104) 


Then multiplying both sides of Eq. (104) by y 2 (t) and dividing by dt result in the same equation as Eq. (100). 
Thus, V (t) is uniformly asymptotically stable and y(t) is bounded by the same equation as Eq. (101). Therefore, y 
given by Eq. (92) satisfies Eq. (104). 

In this study, the comparison lemma allows the stability of nonlinear adaptive control to be analyzed in a local 
sense using its bounded linear approximation. 


The adaptive law (85) thus now can be expressed by its bounded linear version as 

© T (t ) 4> (x (t ) ) = -B T Pe (t ) <I> T (x (t ) ) T<I> (x ( t))<-yB r Pe{t ) (1 05) 

where y = inf (e[f . )( . +ro) (x(r))r<I>(x(T))dT^ >0Gl for t G + T 0 ), where t 0 = 0, t t = f,-_i + T 0 , and 

i = 1,2, . . . ,n — > 00 . 
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C. Time Delay Margin by Matrix Measure Method 


In the presence of unstructured uncertainty, sufficient robustness must be built into the design of an adaptive controller. 
Time delay margin can be thought of as a measure of stability robustness for adaptive control. If the bound on the 
unstructured uncertainty A ( x,u,z,t ) could be determined, then a time delay margin can be estimated for an input time- 
delay system with an equivalent stability behavior. Thus, instead of studying the stability of the original system (77), 
one seeks to investigate stability robustness of the following equivalent input time-delay system 


x(?) = Ax(?) +B u (t — t d ) + (x(t)) 


(106) 


The following two problem statements are equivalent: 1) for a specified time delay margin, what is the largest 
adaptive gain T that can be used in the adaptive law (85) for a stable adaptation, and 2) for a given adaptive gain T in 
the adaptive law (85), what is the estimate of the time delay margin? 

Consider a special case when the uncertainty is almost linear in structure; i.e., 

<I> (x) = x + g (x) (107) 

where g (x) is a vector function whose magnitudes of higher order derivatives are much smaller than that of the first- 
order derivative. 

Then by Taylor’s series expansion 


T>(x) = x + g(x 0 ) + 


dgix 


dx 


1 d 2 g(x) 
(x-x 0 ) + — 


2 ! dx 2 


(x-x 0 ) 2 + ... 




= x + g(x 0 ) + 


dg (x) 


dx 


(x — xq) + 6 (x 2 ) (108) 


x=x 0 


The error equation corresponding to the time-delay system (106) can be derived by substituting the time-delay 
version of the controller from Eq. (81), thus resulting in 

e (?) = A m x m (, t ) + B m r ( t)-Ax{t)-B K x x (t-t d )+ K r r (t - t d ) - © T (? - t d ) T> (x (? -t d )) + ©* t T> (x (? ) ) (109) 

which upon simplification can be expressed as 

e(t)=Ae (, t ) + BK x e (: t - t d ) + Bu ad (t - t d ) - B&* 4> (x (t )) + BK X [x m (t) - x m (, t - t d )] +BK, [r (?) -r(t- t d )\ (110) 

Using the bounded linear approximation of the adaptive law (85), one gets a piecewise locally bounded linear 
approximation of the adaptive law (85) 


«w(0 = ® (0«»(*(0)+® {t)®(x(t))~-rB'Pe{t)+e' (0«*»(jt(0) 

for t G [ti,tj + To), where to = 0, t, = ?,•_ i + To, and i = 1,2, . . . ,n — > ». 

The second term in the right hand side can be locally approximated by a first -order Taylor’s series as 


( 111 ) 


© T (t ) <*> (x (t)) = © T (ti) <i> (x (?,')) + © T (ti) 4> (x (f/)) At + © T (ti) [4» (x (/))-«» (x (?,■))] + . . . 

= ®j4>(x(f))+®j4>iAt + ... ( 112 ) 


where ©, = © (?,) and ©, = © (?,). 

The piecewise bounded linear local approximation of the adaptive law (85) then becomes 

Uad (t) « -yB T Pe (?) + ®J<i> (x(?)) + ©^ i>,A? (113) 

Furthermore, the term <t (x (?)) can be approximated by a first -order Taylor’s series as 


<&(x(?)) = 




dx 


m = 


dx 


(ti)' 


dH > (x(? , '))-x(? I - ) [x(?)-x (?,-)]■ 


dx 2 


d<t>(x(tj )) 
dx 


[x (?)-x (?«)]■ 


= $,■ [x,n (?) - e (?)] + 4>,- Xj [x, n (?) - e (?) - Xj\ + . . . (114) 
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where <fr = <9<t> (*(?,)) jdx and <f> ; = <9 2 <f> (*(?,)) / dx 2 . 

Differentiating the error equation (110) yields 

e (?) = Ae (?) + BK x e (? - t d ) - jBB t Pe (t-t d ) + B®) [x m (? -t d ) - e(t -t d )] 

+ B®J 4>, ii [x m ( t -t d )-e(t- t d ) - x t } + B@J 4>,A? - B ®* T < !>, [x m {t)-e(t)\ 

- // 0 T <f> ; ii [ x m (?) - e (?) - Xj\ + BK X [x m (?) - x m ( t - t d )\ +BK r [r (?) -r{t- t d )\ (115) 

for ? £ [?/,?,- + To), where ? o = 0, ?,- = ?,•_ i + To, and i = 1 , 2 , . . . ,n 

Thus, the locally bounded linear approximation of the error equation can be expressed as 







where 


,4 -I «(-)■ 1 <b' B ® 1 (pjx, 

I 0 


e(?) 
e (?) 


bk x — b®j <i>; 

— yBB T P — B®J <$>"xi 

2? 

1 


d\ (t)+d2(t-t d )+d3 

0 

0 

l 


0 


d\ (?) = ~B@* 0 ( i m (?) — fi©* T <f> ; XjX m (?) + BK x x m (?) + BK r r (?) 
di (? - t d ) = //©, (? - ? d ) + «©/ (?-?</)- -S/Tri'm (? - t d ) - BK, r (? - ? d ) 

c/3 = —B®J <f>, XiXj + B®J <j>,A? + B® iT (\x x,Xi 


(116) 

(117) 

(118) 
(119) 


are treated as time varying disturbances and therefore do not affect the closed-loop adaptive system stability. 

Equation (116) shows that the stability of the locally bounded linear approximation depends not only on the 
adaptive gain T but also on the trajectory of the state vector x(t) together captured in the parameter 7 . Furthermore, it 
also depends on the initial values of the state vector v, and its time derivative i,-, the weight ©,, and the derivatives <t> ; - 
and <t>" , as well as the unknown parameter <f>*. Using the matrix measure method as defined previously, the time delay 
margin of the adaptive system can be estimated as 


(Oi<n(-jCi) + IIAII 


( 120 ) 


where 


t di < — cos 

COi 


A i(Ci)+n(jDj 

II All 


C,; = 


4 4 B® 1 <b' BQ 1 O' x, 

/ 0 


( 121 ) 


( 122 ) 


Dj = 


BK X 4 B®] O' yBB r P + B®J<p'-Xi 

0 0 


(123) 


for ? £ [?,-,?,■ + To), where ? 0 = 0, ?,- = ?,_i + To, and i = 1,2, . . . ,n 

It is noted that the computation of the local time delay margin is retrospective in that the estimate is computed for 
a time window for which most recent data have been collected for analysis. 


D. Estimation of Time Delay Margin for Scalar Adaptive Systems 

Consider a scalar system with linear structured uncertainty 

x(?) = ax(t) + b [ 11 (?) + 0*x(?)] (124) 

The reference model is given by 

x m (?) = a m x m (?) + b„,r (?) (125) 

The controller is given by 

u (?) = k x x (?) + k r r (?) — 0 (?)x(?) (126) 

0 (?) = —Fx{t)pbe{t) (127) 
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The input time delay version of the system is 


x(t) = ax(t) +b[u (; t — td) + 0*x(f)] 

From the response of the delay-free scalar system, the parameter y is evaluated as 

p rti+To 


y=— y xf-(x)dT 


The matrices C, and Z), are 


Q = 


A = 


a + b6* 0 

1 0 

-bk x + bdi yb 2 p 
0 0 


(128) 

(129) 

(130) 

(131) 


The following parameters are computed analytically as 

a + bd* + yj (a + bd*) 2 +1 


Ji(Ci) = 


A* {-jQ) = 
FC/A) = 


2 

1 
2 

yb 2 p 


|| Ail = \bk x -b6i) 2 + y 2 b 4 p 2 

The cross-over frequency and time delay margin are then estimated as 

(Oj < Ji (- jCt ) + 1 1 A|| = \ + \J {bk x - bdi) 2 + y 2 b 4 p 2 

2 ^a + bd* + \J (a + bO*) 2 + 1 + yb 2 p 


t d .< lcos-^( Q) „ + ^ aA 
' (Oi Ai 


1 + 2a/ (bk x — bdi) + y 2 b 4 p 2 


(132) 

(133) 

(134) 

(135) 

(136) 

(137) 


2J{bk x -bei)- + y 2 b*p 2 


The “exact” values of the cross-over frequency and time delay margin for the locally bounded linear approximation 
of the error equation can be determined as follows: 


det (jo* - Q + D i e ~ i<0 ‘ td ' ) = 0 

This results in two equations 

- co 2 - (bk x -bdi) (Oi sin a>,t di + yb 2 p cos (Oit di = 0 

— (a + bd*) (Oi — ( bk x — bdi) (Oi cos (Oit dj — yb 2 p sin a )jt dj = 0 
The cross-over frequency equation is obtained as 


co 


(a + bd*) 2 — (bk x ~ bdi) 2 co 2 - fb A p 2 = 0 


(138) 

(139) 

(140) 

(141) 


The solution gives 


(Oj = 


(bk x - bdi) 2 - (a + bd*) 2 + ^ 

1 4y 2 b 4 p 2 + 

+ 

s- 

CD 

* 

to 

1 

1 

<3- 

cp 

to_ 

2 


(142) 
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( 143 ) 


_ 1 -i Yb 2 P (D? - ( a + b9*) ( bk x - bOj) (0? 

tj. COS r\ 

' y 2 b 4 p 2 + ( bk x — bOj) 2 (Of 

The adaptive gain for which the adaptive system is stable, independent of time delay, can be estimated as 

/ -a-be*+\/(-a-bG*) 2 + l 

1 1 A II = V(-^ + ^«) +7W < -/£(C«) = \ (144) 


Tq \ / 7 + (a + bd*) 


l + y/l + (a + bO*Y 


- (bk x -bdi) 


r< 


b 2 p f t ' +T x 2 (x)dx 


(145) 


provided x(t) ^ 0 for all t £ [ti,ti + 7o). 

Metrics-driven adaptive control is a concept whereby the adaptive gain T can be adjusted for each time window 
in order to achieve a desired or specified time delay margin. This type of adaptation can allow a trade-off between 
transient performance and stability robustness. Let T £ / be a desired time delay margin, then the metrics-driven adaptive 
gain can be computed from 




1 -i /i(Q-i) + /i(jA-i(^)) 

UA-tWII 



(146) 


This is a nonlinear equation which can be solved for the metrics-driven adaptive gain T, for a current time window 
t £ [ti,ti + 7o) using the information from a previous time window t £ [f,_ i , i + 7o). It is noted that the adaptive gain 
r is inversely proportional to the mean-square value of the system state j- J^' j 1 x 2 (t)c/t, which has a notion of the 
system state “energy”. Thus, stability of the adaptive system is dependent on the product of the adaptive gain and 
the system state mean-square value. For metrics-driven adaptive control to achieve a desired time delay margin, this 
product needs to be kept at a desired value. Therefore, if the system state mean-square value is high, then the adaptive 
gain needs to be reduced, and vice versa. 

Example: Given a = 1, b = 1, 0* = 0.1, a m = — 1, b m = 1, p = 1, 9 (0) = 0, r(t) = sin (/). The control gains are 
computed to be k x = — 2 and k r = 1. 

For r = 1, the time histories of the state x(t), control it (f ), and the weight 0 (t) for the delay-free system are 
plotted in Fig. 1. It is noted that the weight 9 (t) converges to the correct value after about t = 40 sec. 



Fig. 1 - Time Histories of Delay-Free Adaptive Control System 
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Select a time window 7o = 1 sec. For the first window 0 < t < 7o, from the response of the delay-free scalar 
system, the parameter y is evaluated numerically at t = To as 

y = — f ° Fx 2 (t) dx = 0.0266 
To Jo 

and 

At (Ci) = 1.293 
P(-;'Ci) = i 
Ji{jDy) = 0.0133 

||£>i || =2.000 

The bounds on the cross-over frequency and time delay margin are calculated to be 

oj\ < 2.500 rad/sec 
td j < 0.344 sec 

For comparison, the “exact” values of the cross-over frequency and time delay margin for the locally bounded 
linear approximation of the error equation are computed to be 

ft)] = 1.670 rad/sec 

t ( [ ] — 0.591 sec 

and for the non-adaptive LTI system for which 9 (t ) =0 for all t are 

ft) = 1.670 rad/sec 
td = 0.592 sec 

The “exact” results are almost the same since the estimates using the locally bounded linear approximation of the 
error equation are for the first time window for whixh 0/ = 0. 

The numerical evidence of the time delay margin is estimated to be t*, ~ 0.407 sec. Thus, the estimated local 
time delay margin for the first time window is in a reasonable agreement with the numerical evidence. On the other 
hand, the “exact” value of the time delay margin for the locally bounded linear approximation of the error equation 
over-estimates the time delay margin of the adaptive system. The process is then repeated for the next time window 
and so on. 

Figure 2 is a plot of the variation of the local time delay margin estimates within the time interval for three different 
sizes of time window; 7o = 1 sec, Tq = 5 sec, and 7 q = 1 0 sec. It is noted that as the window size increases, the variation 
in the local time delay margin estimates decreases. The longer time window allows a more uniform average value of 
the parameter y to be computed, thus reducing the local variation of time delay margin estimate from one time window 
to another. It appears that the mean value of the computed time delay margins is relatively insensitive to the window 
size. In fact, the mean estimate of the time delay margin for the three time window sizes 7o = 1 sec, 7o = 5 sec, and 
To = 10 sec are 0.319 sec, 0.319 sec, and 0.318 sec, respectively. 

Figure 3 is a plot of the mean value of the time delay margin estimates as a function of the unknown parameter 
— 1 < 9* < 1 for To = 1 sec. Generally, 9* is not known, so in a verification setting, time delay margin should be 
computed over all possible parameter variations within their physical bounds. As can be seen in Fig. 3, for 9* < 0, 
the time delay margin estimate of the adaptive system is greater than that when 9* > 0. This is consistent with the 
observation that for 9* < 0, the open-loop system is more stable than that when 9* > 0. Also plotted is the time delay 
numerical evidence from simulations. The numerical evidence is larger than the mean value of the time delay margin 
estimates as computed from the bounded linear stability analysis method by about 20 to 30 percent. Nonetheless, 
the estimation of the time delay margin by the bounded linear stability analysis method is quite reasonable and, more 
importantly, is not too overly conservative. This is important from a practical perspective since any analytical tool for 
estimating time delay margin for an adaptive system must be realistic with reliability, good accuracy, and sufficient 
conservatism. 

Figure 4 is a plot of the mean value of the time delay margin estimates as a function of the adaptive gain 1 < 
r < 100 for 7o = 1 sec. It can be seen that as the adaptive gain T increases, the time delay margin of the adaptive 
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system decreases. This is a well-known phenomenon in the conventional model reference adaptive control. Thus, 
there exists a trade-off between transient performance and stability robustness. Increasing the adaptive gain F gives 
better transient performance but at the expense of robustness. Metrics-driven adaptive control could provide a way to 
maintain consistent time delay margin throughout adaptation in exchange for lower transient performance. 



Fig. 2 - Time Delay Margin Estimates 



Fig. 3 - Time Delay Margin Variation with Unknown Parameter d* 
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0.4 


0.35 - 



Fig. 4 - Time Delay Margin Variation with Adaptive Gain 

V. Conclusions 

This paper presents a new method for estimating time delay margin for model reference adaptive control. The 
bounded linear stability analysis method provides a locally bounded linear approximation of the conventional model 
reference adaptive law. In effect, the adaptive law is transformed into a locally bounded linear approximation within 
small time windows for which local time delay margins are to be estimated. A matrix measure approach provides 
a simple analytical method for computing an upper bound of time delay margin for a linear system is introduced. 
This method is shown to provide a good estimate of time delay margin for a linear system without incurring too 
much conservatism. To analyze the time delay margin, the adaptive system is formulated as an input-time-delay error 
equation. Using this method for the locally bounded linear approximation of the input-time-delay error equation, time 
delay margin for a model reference adaptive control can be estimated. A special case for a scalar model reference 
adaptive control system is studied. The method was able to provide a reasonable, yet not too conservative estimate 
of the time delay margin for the scalar adaptive system. The effect of the time window size was examined. As 
the time window size increases, the variation in the estimates of the local time delay margins for each time window 
decreases. It is found that the mean value of the local time delay margin estimates over the time interval is quite 
insensitive to the time window size. Thus, the mean value of the local time delay margin estimates could be considered 
as a representative time delay margin estimate for the entire time interval. The method also correctly predicts a 
typical behavior in model reference adaptive control whereby the time delay margin decreases with an increase in the 
adaptive gain. Future work will relax the restriction of almost linear uncertainty and also extend this method to robust 
modification schemes in adaptive control. 
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